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Abstract 

The contributions to the muon anomalous magnetic moment from hadronic 
vacuum polarization and from hadronic light-by-light scattering are reexamined 
within the frame work of chiral perturbation theory; the 1/A^c-expansion; and 
the extended Nambu Jona-Lasinio model of low-energy QCD. 
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The anomalous magnetic moment of the muon is known to an accuracy of less than 
lOppm. The latest result from the combined and measurements performed at CERN 
is [1] 

= ^(^m - 2) = 11 659 230(85) x 10-'°. (1) 
The standard model prediction quoted in the most recent review article on the subject [2] is 

^SM ^ ^-^ ggg 192.8(17.6) X 10-'°. (2) 

This includes the pure QED contributions from leptons to four-loops in perturbation theory, 
as well as the lowest order weak interaction contribution 

^weak ^ 19.5(1.0) X 10-'°; (3) 

and the contributions from hadronic vacuum polarization (Figs, la and lb), and from the 
hadronic light- by- light scattering (Fig. 2). These hadronic contributions amount to a total 

= 702.7(17.5) x 10-'°. (4) 

The agreement between the standard model prediction and the present experimental value is 
rather good. 

A new high-precision experiment to measure with an expected error of ±4 x 10"'*^ is 
underway at BNL. To confront usefully theory with the new experiment will require a better 
determination of the hadronic contributions. A substantial improvement in the accuracy of 
the contribution from the hadronic vacuum polarization is possible. It requires, essentially, a 
more precise determination of the low energy tail of the total e+e- — > hadrons cross-section. 
The VEPP-2M facility currently in operation, and DA$NE in the near future, can attain that 
precision. By contrast, the contribution to from the hadronic light-by-light scattering, 
which we shall denote a^ijiad 2), cannot be expressed as a convolution of experimentally 
accessible observables. Here, we have to resort to theoretical wisdom The result quoted in 
[2], based on the estimates made by Kinoshita, Nizic and Okamoto (KNO) [3], is 

a^{had 2) = 4.9(.5) x 10-'° . (5) 

Doubts on the reliability of the accuracy of the hadronic light-by-light scattering contri- 
bution above have been recently raised in refs. [4] and [5] . As emphasized by Martin Einhorn 



[5] , the issue is important for the interpretation of the new experiment and crucial for drawing 
inferences about potential new physics. The purpose of this note is to analyze the hadronic 
contributions to a^^ from the point of view of low energy QCD and to show that, within this 
framework, it might be possible to improve the determination of af^{had 2) to a sufficiently 
safe accuracy for a useful interpretation of the new BNL-experiment. 

1. I suggest we look first at the well known contribution from hadronic vacuum polarization 
in Fig. la which we shall denote a^{had la). It is usually expressed in the form of an integral 
representation (e^ = 47ra, h = c = 1) 

aJhad la) = - T -K{t/ml)e^-ImU^^\t) , (6) 
n Jq t ^ TT 

which is a convolution of the hadronic spectral function —Imll^^^t), related to the total 

TT 

e+e" —>■ hadrons cross-section a{t) {me 0) by 



a(0 = ^e2i/mn(^)(0, (7) 

t TT 



2. ; , X^{l-X) 



with the QED function 

Kit/mi) = / dx ^ — 7/ o , (8) 
^ Jo a;2 + (1 - x)t/ml ^ ' 

which is positive and monotonically decreasing from t = Otoi = oo. I propose we try 
to calculate a^{had la), without using the experimentally known cross-section cr(t). This 
calculation will then serve us as a test of how well we understand - theoretically - a^ihad la). 
(If here we fail to reproduce the phenomenological determination, there is little hope to expect 
that we can make a reliable estimate of a^{had 2).) 

For our purposes, it is also convenient to express a^^had la) using the integral represen- 
tation [6] 



au{had la) = — dx{l — x) 
TT Jo 



2-niH) I _2 



(9) 



where {Q ) denotes the renormalized (H^ '^(0) = 0) hadronic photon self-energy. The 
link between the two integral representations in (6) and (9) is provided by the once subtracted 
dispersion relation (Q^ = —q^ > in our metric +, — , — , — ) 

nk'"(Q^)= Tyr^iwiW. (10) 

Jo t t - TT 



The bulk of the integral in (9) is governed by the low energy behaviour of nSf^(Q2). The 
typical momemtum of the virtual photon in Fig. la is ~ m^. These are momenta values 
much smaller than the characteristic scale of the spontaneous chiral symmetry breaking 
in QCD with three light flavours (A^^ ~ IGeV). Therefore, the appropriate way to look at 
this problem is within the framework of the low energy effective field theory of QCD i.e., 
chiral perturbation theory (xPT). 

Since 11^ '^(0) = 0, the leading contribution to a^{had la) from tree level couplings in 
xPT comes from the 0{p^) effective Lagrangian ; and more precisely from the term propor- 
tional to Ffj^yuF^^ . The coupling constant of this term is related to the slope at the origin 
of the hadronic photon self-energy. The relevant effective Lagrangian we are concerned with 
here is 

£e// = -\ [F^uF^^^ - Pie^a^F^^^AF^, + ...} (11) 



where 



(12) 

Q2=o 



In QCD, and in the large (number of colours) approximation [7], the constant Pi adquires 
a leading contribution of 0{Nc). We can then conclude that to leading order in the 1/Nc- 
expansion of QCD and in the chiral limit, to lowest non-trivial order in xPT, the calculation 
of a^{had la) can be reduced to the determination of only one coupling constant of the 0{p^) 
chiral effective Lagrangian of QCD, with the result [8] 

a^ihad la) ~ (-) — . (13) 

What do we know about the constant Pi? As already shown. Pi is an 0{p^) coupling and, 
unfortunately, our present knowledge in xPT stops at 0{p^). The coupling constants of the 
O(p^) effective chiral Lagrangian of QCD have been determined phenomenologically in a series 
of beautiful papers by Gasser and Leutwyler [9]. It has subsequently been shown [10, 11] that 
these constants are practically saturated by the lowest resonance exchanges and particularly 
by vector-exchange whenever vectors can contribute. We can adopt this phenomenological 
observation to obtain an estimate of Pi . Using the most natural description of vector mesons 
in terms of vector fields as discussed e.g. in ref. [12], one obtains the result 
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(14) 
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where fy denotes the couphng of vector mesons to the electromagnetic field (/y ~ 0.2) and 
My is the vector meson mass {My — 770MeV). The resulting estimate of an{had la) in eq. 
(13) is then 

a^[had la, eqs. (13, 14)] ~ 7.1 x 10"^ , (15) 

to be compared to the latest phenomenological determinations [3, 13, 14] 

{ 7.068(0.059) (0.164) x 10"^ , (a) 
6.84(0.11) X 10"^ , (6) (16) 

7.100(0.105)(0.49) X 10"^ (c) 
where in (16a) and (16c) the first error is statistical and the second systematic. The agreement 
of our xPT estimate with the full phenomenological determination is embarrassingly good. 

At this point some readers will object that the estimate in (15) uses, after all, phenomeno- 
logical input - the ratio fy/My. I want to discuss next, as a possible alternative, how well the 
extended Nambu Jona-Lasinio model (ENJL-model) of low energy QCD, which has recently 
been developped in refs. [15] and [16], does in calculating a^^ijiad la). The ENJL-model can be 
viewed as an approximation of large- A^c QCD, where the only new interaction terms retained 
after integration of the high frequency modes of the quark and gluon fields down to a scale 
at which spontaneous chiral symmetry breaking occurs, are those which can be cast in the 
form of four-fermion operators. The parameters of the model are then A-^ and the two coupling 
constants Gs and Gy of the two-possible (scalar-pseudoscalar) and (vector-axial) four-fermion 
couplings. These two couplings can be traded for the mass Mq of the constituent chiral quark, 
which appears as a non-trivial solution to the gap equation involving Gs', and the efi^ective 
axial coupling qa of the constituent chiral quarks to the pseudoscalar Goldstone bosons [17] 

9a = ^ , (17) 

l + 4G„^r(0,^) 

X X 

where r(0,£) denotes the incomplete gamma function 

T{n,e) = / —e-'z'' . (18) 

Je Z 

In ref. [16] the QCD vector two-point function, within the approximation provided by 
the ENJL-model, has been calculated to all orders in powers of momenta at the leading order 
in the l/Ag-expansion with the result 

nw = "v'W') (19) 

= X 



where 



nv(Q) = Y^8y^ dyy{l-y)T\0, 1 . (20) 



In terms of IIy\ the hadronic photon self-energy is given by the expression (SQ^ — 2/3,i = 
u, d, s) 

n5f^(Q^) = ^(4^(g^)-4^(0)) ; (21) 
and the 0{p^) couphng constant Pi defined in eqs. (11) and (12) can then be easily obtained: 



P.(ENJL)= ^=214 
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(22) 



In the limit Gy ^ (i.e., qa = 1) and >> Mq we recover the result of the constituent 
chiral quark model [18], [19]. Using the input parameter values 

Mq = 265MeV, = 1165MeV, qa = 0.61 (23) 

corresponding to the fit 1 in ref. [15] one gets 

[had la ; eqs. (13, 22)] = 6.7 x 10^^ (24) 

while, using the same Mq = 265MeV value, the constituent quark model prediction comes 
out to be rather low (3.9 x 10~^). We conclude that the EN JL- model does well in reproducing 
the semiphenomenological estimate of eq. (15). 

The deep reason why I am discussing the EN JL- model is that, as I shall next show, it 
offers the possibility to go beyond the leading 0{p^) contribution of xPT. When, later, we 
discuss the hadronic light-by-light contribution, this will appear to be a crucial issue. 

In trying to estimate the 0{p^) correction from the hadronic photon self-energy to the 
integral in eq. (9) there appears an interesting difficulty. Once we factorize the second 
derivative at the origin of the photon self-energy in eq. (9), the remaining integral over the 
a;- variable diverges logarithmically in the ultraviolet region (x ^ 1). In the language of xPT 
this means that a genuine local countcrtcrm of the type ifja'^'^tfjF^i, has to be introduced, with 
a constant which after renormalization will absorb the C/y-divergence of the previous loop 
integral. However, the finite part of this new constant is not fixed by arguments of symmetry 
and depends on the details of the dynamics of the underlying theory. It is precisely in this 
sense that the ENJL-model helps, without appealing to extra phenomenological input. If we 



insert the full expression for the hadronic photon self-energy given by eqs.(21, 19, 20) in the 
integrand of eq.(9), the integral over the x- variable is convergent. Strictly speaking, the Q^- 
dependence of the hadronic photon self-energy predicted by the ENJL-model is expected to be 
adequate up to values smaller than A^. Beyond that we enter a regime where a description 
in terms of perturbative QCD is certainly more appropriate. We have here another observable 
which, like the 7r+ — 7r° electromagnetic mass difference recently discussed in ref. [16] (see also 
references therein), offers the possibility to test the quality of the matching between long- 
distance behaviour and short-distance behaviour. I don't want to discuss this in detail now, 
because it would take us too far away from the main purpose of this note. Here, I shall only 
give the result I get for the total a^{had la) using the ENJL-model expression of n^'*((5^) 
for the evaluation of the long-distance part contribution; and using perturbative QCD for 
the calculation of the short-distance part contribution. The matching between short-and 
long-distances is defined by the optimal choice of an x in the integral (9) which minimizes 

the variation of the total long ( / ) - plus short ( / ) - distance contribution. With the same 

Jq Jx 
input parameter values as in eq.(23), I find 

a^{had la, ENJL ) = 6.7 x 10"^ . (25) 

This is the result corresponding to the optimal choice x ~ 0.9. The stability around this value 
is very good. (Perhaps I should clarify that, in terms of the t-variable in the hadronic spectral 
function representation of eqs.(6, 7, 8), the value x ~ 0.9 corresponds to an "equivalent" 
t ~ 2GeV^; which is very reasonable). The fact that the two results in (24) and (25) are 
exactly the same is a numerical coincidence. If instead of using the values of the input 
parameters in (23), I use the values corresponding to fit 2 in [15] (i.e., Mq = 263MeF, A^^ = 
1048MeF, QA = 0.62) the result is a^(/iad la, ENJL fit 2) = 6.6 x 10"^ 

So far, I have not discussed the explicit contributions to a^(/ia(i la) from tt+tt" (and 
K'^K~) vacuum polarization insertions in Fig. la. In the framework of xPT, these contribu- 
tions appear as chiral loops induced by the lowest order effective Lagrangian 

^eff = \fltTD^UD^U^ + ... , (26) 

where C/ is a unitary 3x3 matrix with det U =1 which collects the pseudoscalar fields and 
DfJJ the covariant derivative 

Df^U = d^U - i{v^ + a^)U + iU{v^ - a^) (27) 



with Vn{afj) external vector (axial) field sources. In our case : = eQA^ and = 0, where 
Q = 1/3 diag (2, —1, —1) and is the electromagnetic gauge-field. The Lagrangian in (26) 
has an efi'ective coupling 



-ie A^\^n^ d^n- d^K- j , (28) 
which leads to the lowest order (i.e., 0{p'^)) hadronic spectral function in xPT : 

3 /2 

^M^foU) = ^l(l-M^ e{t-4ml)+.^K. (29) 

There are two important features emerging here. First, the fact that, by contrast to the 
previous contributions we have discussed, this spectral function is non- leading in the 1/Nc- 
expansion. (This is in fact a rather general property of chiral loops in QCD.) Second, the 
observation that when ^l^'fn^^^oop inserted in the integral representation in eq.(6), although 
nominally it is a leading O(p^) term in chiral power counting, it contributes a very small 
amount to anihad la). Numerically, I find (7r"'"7r~ + K~^K~) 

a^{had 1; xloops) = (0.71 + 0.07) x 10~^ . (30) 

This small result is essentially due to the empirical fact that the physical threshold for 7r"'"7r~ 
(and K~^K~) production is much bigger than the muon mass, therefore suppressing fur- 
ther the, already 1/A^c-suppressed, O(p^) chiral loop contribution. More quantitatively, for 
t > Am^ the function K{t) in eq.(8) is very well approximated by its assymptotic ^rn^/t 
- behaviour. Therefore, the contribution from chiral loops is practically dominated by the 
slope of their contribution to the photon self-energy. The slope in question is 

Pdxloop) = jt^i (4 ^ I?) ■ '''' 

about one order of magnitude smaller than the leading 0{Nc) slope values in eqs.(14) and 
(22). The result in eq.(30) is the leading chiral loop correction, and as such has to be added 
to our full estimate in the chiral limit of the 0{Nc) contribution in eq.(25). This sum gives 
the best theoretical estimate I can think of at present for a^^{had la) - without using the 
empirical knowledge of the e'^e~ — > hadrons total cross-section - i.e., 

a^{had la, Theory ) = 7.5 x 10"^ . (32) 



It is in the high ball park of the phenomenological determinations in (16a, b, c), but not bad 
at all. We conclude that we understand reasonably well, within the framework of low energy 
QCD, a^ihad la). 

2. It is time to examine now the contribution from the hadronic light-by-light scattering in 
Fig. 2 i.e. a^^{had 2). Again, the typical momenta of the virtual photons in Fig. 2 is ~ m^; 
and the appropriate framework to analyze the problem is the one of x^T. At first sight, the 
leading contribution to an{had 2) from direct tree level couplings appears to come from the 
0{p^) chiral effective Lagrangian with four electromagnetic field strength tensors i.e., the 
Lagrangian (SQf = 2/9, i = u, d, s) 



This is the QCD-analogue of the famous Euler-Heisenberg effective Lagrangian of QED. [20]. 
(For an excellent book on the subject see ref. [21].) With the normalization I have chosen, 
the coupling constants hi and h2 are dimensionless. In the constituent chiral quark model 
[18,19,15] hi = h2 = 1. In QCD these effective couplings appear to leading 0{Nc) in the 1/A^c- 
expansion. (This is why I have already pull out the A^c-factor in (33).) The direct insertion in 
the muon vertex diagram of Fig. 2 of the local light-by-light scattering couplings given by the 
effective Lagrangian in (33) leads however to C/V^-divergent integrals. The situation here is 
rather similar to the one we already encountered in the case of hadronic vacuum polarization 
when one tries to insert the local 0{p^) photon self-energy coupling. Again, what we find 
now is that we need more information about the underlying QCD dynamics than just the 
two local couplings hi and /i2- In view of its phenomenological successes it seems reasonable 
to expect that the ENJL-model we discussed earlier may provide the required help. More 
on that later. First, let us examine other possible contributions to a^{had 2) within the 
framework of xPT. 

There is of course the contribution from the hadronic light-by-light scattering am- 
plitude induced by the effective lowest order coupling in (26). Like in the case of the 
hadronic vacuum polarization, it gives rise to an 0{p^) contribution, which is also non- 
leading in the 1/A^c-expansion. The chiral loop calculation of the light-by-light scattering 
tensor H^^'^^'^ {k,qi,q2Tq'i) leads to a rather complicated expression. In order to get some 




167r2 360Mj9^ 
+7/i2(£'^"^"F^.(^ 




(33) 



feeling, we can however consider the hmit where all the photon momenta are negligible with 
respect to the physical pion mass. This limit is described by the Euler-Heisenberg effective 
Lagrangian of scalar-QED [21] i.e., 

£e//(scalar-QED) = -\f,.F^^ + J^^^^^'^^^'^^^aFm. 
+ 1^4'' + {e'^'^F.^F.^f] + ... (34) 

As already discussed earlier, the second term in the r.h.s. above governs the slope of the 
photon self-energy induced by the lowest order chiral loop (see eq.(31)). The third term is 
the resulting one for the effective light-by-light scattering. From the comparison between 
(33) and (34) one can see that for values /ii ~ /i2 ~ 0(1), the order of magnitude of the 
direct 0{Nc) couplings are about the same size (or smaller for the [FFY - term) than those 
induced by the pion chiral loop. The moral here is that the contribution to a^{had 2) from 
chiral loops, although 1 /A^^-suppressed, may well give rise to a sizeable contribution. 

There is yet another contribution to hadronic light-by-light scattering which in xPT 
appears as a tree level contribution to 0{p^) in the chiral expansion; but contrary to the 
0{p^) contribution which we have just discussed, it is leading in the l/N^ - expansion. It is 
the one induced by the Adler [22] , Bell and Jackiw [23] anomaly which in xPT is incorporated 
in the so called Wess-Zumino-Witten term (WZW) [24,25]. When restricted to two external 
photons and one pseudoscalar field, the WZW-term has the simple form 

Ubj = (^n\x) + -^vix)^ . (35) 

The 7r'^(?7) propagator between two effective ABJ-verticcs brings down the resulting 0{p'^) x 
0{p'^) light-by-light scattering amplitude to O(p^); and because of the factor 1//^ the overall 
contribution is 0{Nc), and therefore leading in the 1/Ac-expansion. Obviously, this contri- 
bution cannot be neglected; and in fact, the simple power counting arguments above tell us 
that it could be the dominant contribution. 

The previous analysis of a^{had 2) within the framework of xPT has served to identify 
three distinct sources of a priori important contributions. We want now to discuss their 
possible evaluation. 

The contribution to a^{had 2) from the lowest order chiral pion loop can be calculated 
unambiguously. It has in fact been evaluated numerically by KNO using two different meth- 



ods, with the average result 



ai^{had 2, x pion loop ) = -5.12(0.35) x 10"^° . (36) 

To this result one should add the corresponding contribution from the chiral kaon loop; 
presumably rather small. 

A straightforward calculation of the contribution to a^{had 2) from the 0{p^) light- 
by- light scattering amplitude resulting from the 7r^(r7)-exchange between two AB J- vertices 
leads to t/y-divergent integrals. Once again, this is an indication that more about the 
underlying QCD dynamics than just the local 7r° (77)77 0{p'^) coupling is needed. KNO, 
in their evaluation of the "7r°-exchange" contribution, have solved the problem by using a 
vector- meson-dominance (VMD) inspired regularization of the three photon- propagators; i.e., 
they use the prescription = —qf i = 1, 2, 3) : 

-1-11 

q!^q!^qI + m^ ■ ^^^^ 

It turns out that this VMD - prescription is in fact very close to the result one obtains in the 
ENJL-model [16], after summing all orders in powers of momenta at the leading order in the 
1 /Nc - expansion, which leads to the effective replacement : 

-1,-1 1 _ -Mm) (38) 



with U^y\Q^) as given in eq.(20), and Mv{Q'^ = 0) = My. Numerically, the result found by 
KNO from the "7r°-exchange" term, and using the VMD-prescription in (37), is 

a^(7r°) = 6.5(0.6)10-^° . (39) 

The corresponding contribution from //-exchange which should be added to (39) is likely to 
be negligible. 

An improved calculation of the "tt*^ -exchange" contribution within the framework of the 
ENJL-model is possible. Schematically, it corresponds to the replacement of the yr'^-exchange 
tree level diagram in Fig. 3a by the one in Fig. 3b. The round dots (o) in this figure 
represent four-fermion vector like couplings; the square dots (□) four-fermion axial-vector 
like couplings. The horizontal string of bubbles in Fig. 3b has already been calculated (see 



section 3.2 of ref. [16]). It contains in particular the pion (eta) - pole(s), but also axial- vector 
contributions. The three-point functions at the vertices of Fig. 3b; i.e., the links xQc ^^d 
□ Oo 5 have non-trivial Qj - dependence. In KNO's calculation, they have been approximated 
by the corresponding local AB J- vertices in (35). 

Finally we turn our attention back to the leading 0{Nc) contribution we first discussed, 
which in xPT starts at 0{p^). This is an independent contribution which has to be evaluated 
separately; and, contrary to what has been done so far in the literature, it should be added to 
the two other contributions previously discussed. KNO's numerical evaluation of what they 
call the "quark- loop" contribution provides in my opinion a good first approximation to the 
required calculation. The constituent quark-loop graphs of refs. [3] and [26] correspond to 
the mean field approximation of the ENJL-model of refs. [15,16]. In the case of hadronic 
vacuum polarization, this is the approximation where Hy^ — > 11^^ (see eqs. (19) and (20)). 
The result found by KNO for the "quark- loop" contribution is 

ttf^ihad 2, quark-loop ) = 6.0(0.4) x 10"^° . (40) 

It corresponds to the choice = Md = 300MeF and Ms = 500MeF. (There is also a small 
contribution from the charm quark with a mass value = i-5GeV which has been included 
in (40).) 

Again, it is possible to do an improved calculation of the "quark- loop" type contribu- 
tion within the framework of the ENJL-model. The suggested improvement is schematically 
represented in Fig. 4. In practice it amounts to the effective replacement indicated in (38) 
for each photon-propagator. Since here we are working in the chiral limit the u, d and s 
constituent quark masses should be taken equal. I would expect a sizeable reduction of the 
"quark-loop" estimate in (40) by perhaps as much as a factor of two. 

The main conclusion which emerges from our analyses is that the three contributions in 
(36), (39), and (40) correspond to well identified different sources in low-energy QCD. The 
three contributions should be added in order to get an appropriate estimate of a^{had 2). So 
far, only the result from the chiral 7r-loop contribution in (36) can be considered to be fully 
reliable. (In that respect I agree with the analyses of ref. [5] .) Both the "7r°-exchange" and the 
"quark-loop" results in (39) and (40) should be considered as first approximation estimates; 
specially the "quark- loop" result (40). These estimates can however be substantially refined if 
one is ready to accept the ENJL-model as a good approximation to low-energy QCD at large- 



Nc- We see no reason why the ENJL-model should fail here, and the results we have obtained 
in the case of the hadronic vacuum polarization are very encouraging. The improvement 
can be done following the suggestions I have sketched above. The implementation of this 
improvement is quite a formidable task, but it should be relatively easy to incorporate in the 
already existing computer programs developped by Kinoshita et al. A parallel calculation of 
the 



77 ~^ 11 cross-section at low energies in the ENJL-model, and its comparison to future 
experimental results of DA$NE, could help to check the validity of the model. 
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FIGURE CAPTIONS : 

Fig. 1 Hadronic vacuum polarization to the second order muon — 2 (a) ; and to the fourth- 
order muon g — 2 (b). 

Fig. 2 Hadronic hght-by-hght scattering contribution to the muon g — 2. 

Fig. 3 Hadronic hght-by-hght scattering contribution induced by the combination of two 
three-point functions. Graph (a) represents the 7r°-exchange contribution. Graph 
(b) the corresponding generahzed contribution within the ENJL-model. The square 
dots (□) represent four-fermion axial- vector like couplings ; the round dots (o) four- 
fermion vector like couplings. 

Fig. 4 Hadronic light-by-light scattering contribution induced by the irreducible four-point 
function. Graph (a) represents the mean field approximation. Graph (b) the full 
ENJL-model contribution to all orders in momenta for the off-shell photons. 
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